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Annual Summary for NASA Interchange NCA2-543 


Three technical reports were submitted this year, including this one. 
They are: 

1. Electromagnetic characterization of conformal antennas (Univ. of 
Michigan Report #027723-4-T). 

2. Scattering and radiation analysis of three dimensional cavity arrays 
via a hybrid finite element method (Univ. of Michigan Report 
#027723-5-T). 

3. A finite element-boundary integral formulation for conformal arrays 
on a circular cylinder (Univ. of Michigan Report #027723-6-T). 

The first of these reports described several extensions and further 
developments associated with the microstrip patch array code developed last 
year. Specifically, modifications were introduced into our original finite 
element formulation to permit simulations of embedded resistive cards, 
lumped loads and impedance surfaces. The new implementations were 
validated, and a number of new design schemes were examined for 
controlling the RCS, resonance frequency and array aperture size. For 
example, the RCS gain of rectangular patches was examined as a function of 
lumped load values and locations; dielectric coatings were examined for 
controlling the patch RCS; resistive ribbons were employed and found 
suitable for broadband RCS reduction without excessive compromise in 
gain, dual cavities were shown to provide resonance control and decrease in 
array aperture size; and reactive resistive sheets provided an attractive 
alternative for controlling the patch's resonance frequency. Needless to 
mention, during the course of this year our finite element rectangular patch 
array code was improved in many respects. Its geometry interface was 
expanded to allow for greater adaptability, several new features were added 
as noted above and a variety of reference calculations were generated 
providing potential uses with possible new application of the code. Finally, 
a short code manual was written. 




In addition to the above improvements of our finite element planar 
patch array code, one additional code was completed for infinite planar patch 
arrays, and we began the development of another code suitable for 
cylindrically conformal patch arrays. Specifically, a rather thorough 
analysis and investigation was completed for the scattering and radiation 
analysis of an infinite path array backed by rectangular cavities. The 
pertinent finite element formulation was developed ab initio, and measured 
data were collected for validating the code. By invoking Floquet's theorem, 
the computational domain was restricted to a single element, and this was 
the primary reason for considering the infinite array formulation. We have 
already found that the infinite array approximation is sufficiently accurate 
for modelling large arrays. Thus, instead of using our original patch array 
code which is computationally intensive for large arrays, one can resort to 
the more efficient infinite array code. 

Much of our effort during the second half of this fiscal year was 
devoted to the development of the necessary formulation for the analysis of 
arrays on cylindrical surfaces. The details of this formulation are given in 
the present report following this summary. As of the moment, the entire 
formulation has been documented, and the finite element matrix elements 
were generated using edge-based pie-shell elements conformal to the 
cylindrical aperture. Most of the effort, though, was devoted to the 
discretization of the boundary integral for terminating the mesh. More 
specifically, various asymptotic forms for the cylinder’s Green’s functions 
were examined, and one was proposed for its efficient evaluation without 
compromising accuracy. The implementation of the proposed cylindrical 
array formulation is currently in progress. 
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Abstract 

Conformal antenna arrays offer many cost and weight advantages 
over conventional antenna systems. In the past, antenna designers 
have had to resort to expensive measurements in order to develop a 
conformal array design. This is due to the lack of rigorous mathemat- 
ical models for conformal antenna arrays, and as a result the design 
of conformal arrays is primarily based on planar antenna design con- 
cepts. Recently, we have found the finite element-boundary integral 
method to be very successful in modeling large planar arrays of arbi- 
trary composition in a metallic plane. Herewith we shall extend this 
formulation for conformal arrays on large metallic cylinders. In this we 
develop the mathematical formulation. In particular we discuss the fi- 
nite element equations, the shape elements, and the boundary integral 
evaluation, and it is shown how this formulation can be applied with 
minimal computation and memory requirements. The implementation 
shall be discussed in a later report. 


1 Introduction 

Conformal antenna arrays are attractive for aircraft, spacecraft, and land 
vehicle applications since these systems possess low weight, flexibility, and 
cost advantages over conventional antennas. The majority of previous devel- 
opments in conformal antennas has been conducted experimentally due to a 
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lack of rigorous analysis techniques. Various approximate analysis techniques 
are restricted in many respects, including accuracy and element shape, and 
are based on planar antenna models. 

Recently, we have found that the finite element-boundary integral (FEM- 
BI) method can be successfully employed for the analysis of large planar 
arrays of arbitrary composition [1], The resulting system is sparse due to the 
local nature of the finite element method whereas the boundary integral is 
convolutional, thus ensuring an 0(N ) memory demand for the entire system. 

In this report we will extend the FEM-BI formulation for aperture an- 
tennas conformal to a cylindrical metallic surface. Both the radiation and 
scattering problems will be developed in the context of the FEM-BI method. 
In contrast to the planar aperture array, the implementation of the cylindri- 
cally conformal array requires shell shape elements rather than bricks, and 
the required external Green's function is that of the circular perfectly con- 
ducting cylinder. In its exact form this Green’s function is an infinite series 
which must be evaluated efficiently and must also be put in a convenient 
convolutional form for storage minimization. 

This report presents the FEM-BI formulation, appropriate cylindrical 
shell elements, and the system evaluation strategy which will maintain low 
memory and computational load. The cylindrical elements will be chosen di- 
vergenceless while maintaining excellent geometrical fidelity. These elements 
are derived using the procedure attributed to Whitney [2], Substantial effort 
is also devoted to the development of expressions for the metallic cylinder’s 
dyadic Green’s function which are convenient for computation, and extremely 
inexpensive asymptotic evaluations will be derived. 

2 FEM-BI Formulation 

Consider the configuration illustrated in Fig. 1 where a cavity is situated on 
a circular metallic cylinder. The cavity is recessed in the cylindrical surface 
and its walls shall be assumed to coincide with either constant o or constant 
~ planes. Also, it shall be assumed that radiating elements may reside on 
the surface of the substrate within the cavities. 

In accordance with the FEM-BI formulation discussed in [1], the radiated 
or plane scattered fields can be determined by considering two systems of 
equations. Particularly, the fields interior to the cavity are formulated by 
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the finite element method which results in a sparse system of equations. As 
usual, the determination of the finite element mesh at the aperture requires 
the imposition of an externally supplied condition. For our case this condition 
becomes the boundary integral equation enforced at the cavity aperture. 

We begin the development of the FEM-BI system by first discretizing the 
weighted vector wave equation in the interior of the cavity using cvlindrical 
shell elements which are most appropriate for this geometry. Assuming the 
presence of possible sources and A/' in the cavity the vector wave equation 
is of the form 


V x 


V x E(p. q>, z ) 


pAp-O- z) 

—jk' 0 Z 0 J l (p,Q.z) + V X 


- klt r {p, 0 ,z)E(p.e,z) 


A/ 1 (p, <?>,-) 

Pr(p,<P,z) 


(i; 


In this E(p,o.z ) denotes the total electric field in the usual polar cylindrical 
coordinates. f,.(p.o.z) and p r (p.o,z) are the relative permittivity and per- 
meability of the substrate or material filling the cavity, k 0 is the free space 
wavenumber and Zo denotes the free space intrinsic impedance. To generate 
a system of equations from (1) we apply the method of weighted residuals. 
We obtain 


1, 


V x 


V x E(p, 0 , z) 


iF,(p, 0, z)pdp dcp dz 


fl r (p,0,z) 

~k 2 , f t r(p- d. z) E(p. 0. z) ■][',( p. O. z)p dp do d: 

J ^ i 

L v * 

J ^ O Z a 


O- -) 

P, (P. O.Z) 

J v J‘(p,<?,z) 


■ 11 dp. o, z )p dp do dz 
) • W,(p,0,z)pdpd0dz 


( 2 ) 


where H t (p. 0, z) is a subdomain vector- valued weight function to be specified 
and 1 ; is the fth volume element resulting from a discretization of the cavity. 
Given the sources ( J‘..M ') for each weighting function the right-hand size 
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of (2) is known and shall be referred to as the excitation function defined bv 



.U'(p, o, ;) 

ft T (p, 0 .z) 


j k 0 Z 0 J\p, <p , c) 


W,(p, O. z)pdpdo dz( 3) 


Upon application of a standard vector identity and the divergence theo- 
rem [3], we recognize (2) as the weak form of the wave equation 


k 


V x E{p< d>, ~) • V x \Vi{p, <p,z) 


pdp d<f> dz 


Pr(p , d>, ?) 

■k* / t r (p, 0 ,z)E(p, 0 , z) ■ Wi(p, o , z)pdpd(pdz 

Jv. 


-jk 0 Z 0 j> ii(p, 0 ,z) x H(p, 4 >, z) • Wi{p, 0 , z)dS = /, (4) 


with h(p,<j>,z) indicating the outward pointing normal of the zth element. 
Si is the surface area of that element, and H(p,4>,z) is the total magnetic 
field. It can be shown that the surface integral of (4) vanishes for all those 
elements which do not border the cavity aperture. Furthermore, their non- 
zero contribution is limited to the portion of their surface which coincides 
with the aperture. Thus. S, is a subdivision of the aperture surface, and the 
surface integral can then be more specifically written as 

I p(a.0.z) x H(a,0,z) ■ W t (a,S,z)(ad 0 dz) (5) 

where S a , denotes the /th element of the aperture surface and p(a , d>, c) is the 
unit vector normal to the aperture surface. 

To eliminate H from (4) we introduce the boundary integral equation 

H(a,0,z) = H'(a.0,z) + 

jk 0 Z 0 J p(a, o' , z ) x E(n, 6. z ) ■ G' e2 (u, d>, r; a. 0, z )(ad<j> dz' ) (6) 


which provides an additional relationship between E and H on the aperture. 
In this H'{a.0,z) is the incident magnetic field evaluated on the aperture, 
\ Q = d- is the free-space admittance, and G e2 (a, <j>, z; a, <p', z) is the electric 
dyadic Green’s function of the second kind for a metallic cylinder [4]. This 
Green's function satisfies the radiation condition and the boundary condition 

V x & e i(a,0,z; a, <p,z) = 0. (7) 
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Upon inserting (6) into (4) along with (5). one obtains 
V x E(p, 0 , z) • V x \V,(p, 0, :) 


L 


- p dp d0 dz 


pAp.O.:) 

-kl J t r {p, 0, z)E(p, <t>, z) ■ w,(p, O, ~ )p dp do dz 
+ {k 0 a) 2 1 p{a,0,z)x p{a, 0 , z ) x E(a, <f> , z )• 

*•' Sq i J a 


G e 2 i a ' d. z: a* 0 . z )do dz 


\\\ ( a , 0, z)dddz 


J\ + jk 0 aZ 0 / p(a, o. z) x H l (a, o, z ) • \Y t (a. o . z) 


do dz 


( 8 ) 


which is an equation only in terms of the electric field in the cavity volume 
and on its aperture. It is important to note that S at denotes integration over 
the zth surface element coinciding with the aperture whereas S a indicates 
integration over the entire aperture. 

Following the principles of Galerkin's method for a solution of E appear- 
ing in (8), we expand E in terms of the vector- valued weight functions also 
used for testing, i.e.. 

E(p.o.z) = Y. EjWj(p.o.z). (9) 

j= i 

In this expansion N v is the total number of unknowns or edge fields (interior 
+ aperture edges) and Wj(p, 0,z) are the subdomain vector- valued basis 
functions. By necessity the aperture field takes the form 




where 


E(a , d, z) 


&*(j) 


Y E } 6 a [j)Wj(a,0,z) 

j = i 


1 i f II ^ D apt r tui 


( 10 ) 


= 0 tlst 

Combining (S), (9) and (10) we obtain the FEM-BI system 
v x Wj(p,4>,z) • V X I Vi(p,<j>,z) 


;n) 




p r (p,0,z) 


-pdpdodz 
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~^0 / e r(/>, <5, ~) • I V r ,(/9, 0 , z)pdp d(pd: 

J \ x 

+ {k 0 a) 2 6 a (j) I H , ( a , o. z) ■ p(a,<p,z)x 

Js a , [ 

/ G e2 {a. o. z\ a, <p , - ) x p(a, , c ) • W } ( a , <p , ’ )o?(/> 

^ ^aj 

= fi + i k 0 aZ 0 f p{a,<i),z) x H'{a,4>,z) - W x (a,<p,z)d(i>dz (12) 

J i 

Below we discuss the specifics of the weight/expansion functions. 


dodz 


3 Vector Weight Functions 

To explicitly compute the matrix elements resulting from (12), we must first 
specify the vector- valued weight functions \V t . Traditional node-based shape 
functions associate the system unknowns ( E s ) with the field at a node. In 
contrast, edge-elements have their degrees of freedom associated with the 
field along an edge ol the element. Most importantly, the first order edge- 
based shape functions can be chosen to be divergenceless, thus satisfying an 
inherent characteristic of the unknown field. Comparatively, the first order 
node-based elements are not divergenceless, and hence a penalty function [5] 
must be used to ensure a valid solution of (1). Whitney [2] developed a 
formalism from differential geometry which allows the generation of edge- 
based elements from traditional node-based elements. Although Whitney’s 
procedure does not guarantee divergenceless elements, we will choose a basis 
which is divergenceless. 

\\ hitnev developed a family of p-forms (where p indicates the order of the 
form) from differential geometry which possess characteristics that are use- 
ful to the finite element community. Node-based elements and edge-based 
elements correspond to Whitney 0- and 1-forms, respectively. Traditional 
first-order Lagrange elements (0-forms) provide field continuity which is not 
physical. Namely, these elements have continuity of both tangential and nor- 
mal fields at the element junction whereas physical requirements allow for 
a discontinuous normal component. Edge-elements (1-forms) have only tan- 
gential continuity and are therefore better suited for electromagnetics appli- 
cations. It should also be mentioned that the overspecification of continuity 
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by the O-forms causes spurious solutions of (12) [6]. Although first-order La- 
grange elements are used as the 0-forms in the literature [6], these are mostly 
applicable to elements such as bricks or tetrahedra, and the corresponding 
1-forms are divergenceless by virtue of the first-order 0-forms. However, in 
the case of cylindrical shell elements (see Fig. 2), the 0-form elements cannot 
be of first order if their corresponding 1-forms are expected to be divergence- 
less. To ensure the divergencelessness of the 1-form elements, the 0-form 
cylindrical shell elements are chosen as 


A](p, 0,") — +i 


[fa(P~ Pb)(<f> - 4 >t){z - -() 


A 2 (p, o, r) = - 
A 3 (p, <?, r) = + 
A 4 (p, e>, :) = - 
X 5 (p.O,:) - - 
A 6 (/>, 0. z ) = + 


[pb (p - i°q)(^ - ftrM- - -t) 
p A 

[pb {p- Pa)i0 - <*><)(- - -t) 
P A 

[K (p_z ~ ~ Zll 

P A 

[Pa {_P_ ~ Pb){<t> ~ ~ -6) 

P A 

fpb( pj2 ~ ~ ~ fc ) 

A 


A r(p.o.c) = 

A s(P- o.z) = +i 


[Pb (P ~ Pa){<t> ~ <M(~ ~ ~fc) 
P A 

fp^{ p_2 Pt)(0 - Of )(- ~ ~t) 
A 


(13) 


with A = (p a — pb)(0i — <Pr){~b — z t) and is understood that the support of 
each 0-form element is only within a single shell. The subscripts correspond 
to local node numbers and the element limits (p a , pb- 01, <?> r , -6, ~t) are shown 
in figure 2. An example of a 0-form is shown in figure 3 where \\(p, 0. z) is 
graphed over the three faces of the element. We remark, however, that for 
our application p a and pb are very large and thus for all practical purposes 
the shape elements ( 13) are first order. 

The Whitney 1-forms are generated from the 0-forms by employing the 




relation 


n\j(p,0,z) = \j{p,o<:)V\,(p.d.:) (14) 

where (i,j) refer to the node numbers which define each edge. Denoting the 
normalized shape elements by 


Vij = 


^ ij [Pi Pi ~ ) 
max {\\Wij{p,<j>, s)||} 


From (14) we obtain 


Nuip.e,:) = /w = M» - *)H: - . 

P(<Pr - Ol) 2 {:, - : b ) 2 ; ' 


PDPr (C, - ; b ) p( Qr _ 2(- t _ :b }2 

rise (/>,*, r) = ~ ~ °; } 2 2(C ~ =b)2 p AV(p <P -) = J gq(^-^) 2 (~~--6) 2 

/»( 0 r - 0 /) 2 (-~, - Z b ) 2 ~ — — 


p{<Pr - 0i) 2 {: t - z b )‘ 


-Vh(/3, O.z) = 


- V 23 (/»,<?--) = 


a 2 D> - /^) 2 (~ - 
/’Vt - /M 2 Uf - 

~ da) 2 (- ~ 

P 2 (Pb ~ Pa) 2 (: t ~ 


~i° V 58(/J,0, 


-V(37(P, 0. C) = 


_ />£(/> ~ PbYi = ~ -ft) 2 - 
P 2 (Pb - />„) 2 (r, - r 6 )^° 

_ pKp - ^) 2 (- - -«>) 2 * 


P (Pb - Pa) 2 (=t - -fc)‘ 


N 15 (P,0,Z) = Pb) ! {<:> 0r]2 Z N 26 (p 6 Z) = ^^-Fa) 2 (g>- d> r) 2 . 

^(^6 /0„ ) (<?,• — d>/) 2 ’ ’ p{pb — pn) 2 (0r — 4>t) 2 ~ 

nmp.o.z) = £</>-»)»(*-*)% v ( o = P >[p - P «n* - *,? 

P(Pb Pu ) 2 ( Or O/) 2 p(pb - p u ) 2 {0r - 0l) 2 ' 

Tliese can be put in a more compact form by introducing the definitions 

•V ’,(p.o, S :p.lz) = 

P(<Pr - <A/) 2 (c, - Cft)' 2 

N o (p.0.z:p.lz) = i d ~ P) 2 (P ~ P)H= zl)l 

P 2 {Pb ~ Pa) 2 (=t ~ -ft) 2 

= i'l ~ P)iP ~ P)Ho - i) 1 

p(Pb ~ Pn) 2 (0r ~ 0l) 2 


S 




where d = p a + p b and (p, O, i) are parameters which define each 1-form. 
Comparing (17) with (16), we have 


A^12 ( P, ^ ) A pi Pi *Pi * . * i d*r , ~t ) , 

Ase( Pi ~ ) A P ( p, 0r , ~fe)i 

A'n(d* <f>, = ) = A r 0 (p, <p, pb , •, ~t), 

^23(/9, 0, - ) = Ad(/>. <?> da. •• -(), 


ATj(p, <P , -) = Ap(d> 0. •> -t) 

A r 8-(^ <£, U = A r „(p, <£, •, 0 t , z b ) 

A 58 {p,4>,~) = A r ^(/J, 0, Z\ Pb, Cfc) 

A' 6 7(p, <£, c) = A r „>(p, O, p a , •. C t ) 


Afis(d, <M) = N : (p,0,;\pb,0r,-), A ; 2 6 {p,0,z) = N;(p, 0,Z\ p a ,0r,-) 
N 4 s(p,<t>,z) = N ; {p.0,z; p b ,0 h -), N 37 (p,<f),z) = N : (p,0, z: p a ,0i,-) (18) 

These 1-forms will be used as the vector-valued weight functions to form the 
system (12). Figures 4-6 illustrate some 3-D plots of the weight functions for 
each generic form (17) where we recognize the - p roll-off in the p-component 
as required by the divergenceless condition. 


4 Matrix Elements 

Above, we have derived a FEM-BI formulation for apertures on a circular 
cylinder (12) and developed the appropriate Whitney 1-forms (17) for the 
solution of the system. In this section, we will compute the matrix elements 
associated with the finite element submatrix of the system. The evaluation 
of the elements associated with the boundary integral will be discussed in 
Section 5. 

The overall matrix system resulting from (12) can be symbolically written 
as 


[■M ; is.-,] 

[Co] : [A,] 


(A! = (/,) 


(19) 


The submatrix [.4 U ] is sparse and associated only with the interaction of 
edges interior to the cavity volume and excluding those on the aperture. 
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We shall refer to these edges as belonging to the set J. Whereas the 
aperture edges shall be put in the set A. The sum of the edges in the sets I 
and A include all edges resulting from the discretization of the cavity and its 
aperture excluding those which lie along metal. With this identification, the 
submatrix [Cq] is associated with the interaction between the interior and 
aperture edges, and is also sparse. The last submatrix [D I} } results from the 
boundary integral equation and provides the interactions among the aperture 
fields/edges. It is consequently fully populated, and in the next section we 
shall put it in circulant form to reduce its storage requirement to O(.Y). 

From ( 12 ). the matrix elements can be explicitly written as 


. 4 , 


-k 


V x \\’j(p,Q. z) ■ V x \\’ t {p,<p,z) 


Bi 


- k 


p r (p,o,z) 

-kltr(p* 0- z)W 3 (p. 0, c) • W t (p, o, z) pdpdcpdz {ij G 1 } 
V x lTj(/),o. c) • V x Wi{p, <j>, z) 


p r (p'0,=) 


~ k*t r {p,o,:)\Vj(p , 0.=) ■ if'dp.o, c) 


pdpdo dz 


+{k 0 a) 2 6 a (j) J U ',(a. 0, z ) • 


G e2 (a 

Ci, 


p{a,<p,z)x 


t t , . , / . f 


0, z\ a, 0 , c ) x p(a, 0 , c ) • li'Ja, <j>,z)d<!> dz 
V x Wj(p. 0. z) • V x W,(p,0,z) 


d(pdz > {? G T, j € ^ 4 } 


- i 


pAp.o,:) 


A, 


-k 2 t r (p,0,z)\Vj(p.0.z) ■ iT'dp, o, z) 

-k 

- ko e r ( p. 0, - ) H j ( p, 0 , z ) • Wi {p,4>,z)\pdpd4>dz 


pdpdodz {/ G A. j G X) 


V X O. c) • V X ll'(p, <*>,-) 

Pr{p,0,z) 


+ (k 0 a) 1 6, l (j) / lT'-(fl,o, r) 


p(a. d.z)x 
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a, <j> , 


p(a, o',: ) ■ Wj ( a , 0, 




{ij € A} (20) 


Since w i and ll j are subdomain basis functions whose support is restricted 
to a single element, it is easily recognized that the submatrices [A tJ ] and [C, 7 ] 
are sparse. 

We may now proceed with the evaluation of matrices [A tJ ] and [C tJ ] since 
these do not require a boundary integral. If we assume that the material 
properties are constant within each element, we may define two families of 
integrals which will specify the submatrices. 

The explicit determination of the elements A ZJ and C XJ involves the eval- 
uation of the two integrals. 

Jab ~ j v V x \ a (p.o.z: pi.Oi. ~i) ■ V x N b (p,o.:\p k .O k ,z k )pdpdodz( 21) 
j[l ] = N a (p.O.:: />/.<?/, i/)- N b (p.o.z; p k ,0 k , z k )pdpdodz (22) 


in which (a, b) (E [p. o, ;}, ( l,k ) 6 and N represents one of the generic 

1-forms (17). Carrying out the required vector operations and organizing 
each integral in a separable form we have 


/<’> 


4 2 f Pb dp f*' (O - Ot) 2 (o - Ok) 

J J Pa P J<t>l 

r 

Jp„ P } Jo, 


{Or - Ol )•' 

(o - Ot)(o - o k ) 


(Or ~ Ot) 4 


i0 \„ 

dO / 

JZi 


(r — ~i){: - z k ) 


( ~b — -f) 4 
i.V-il - 


dz 


- : k ) 2 
(=b~=t) A 


d: 


(23) 


! {l) 

pO 


2 Pa 


(<l~Pk ) 2 f Pb {P ~ Pk)(p + Pk) 


L 


L 


(Pb ~ Pa)~ J Pa 
* r [01 - 0) 


dp 


Ol (Or ~ Ol ) 2 


-do 


j: 


i.\2 /_ _ - \2 


(:-m = 


~k J 


( -6 ~t )** 


-d: 


(24) 




-Pa 


(d - Pk) f Pb ip ~ pk){p + pk) 


(Pb - Pa)' 


l 


r<Pr 

J(t>i 


(o - Ot) 2 (o - o k ) 


(Or - Ol) 


^ r ^ 


dp 

it) 


dz 


(25) 


11 




4 ( "- 

Pi?(d- 

Pk ) ( 

Ol) I 

rPb (p ~ 

Pi) 2 (p- 

h)2 do r 

(=-=i)(=- 



( Pb - 

PaY 


3a 

P 3 

dp L 

( z b — Z t 

r d - 

+4^- 
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( 29 ) 
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{ Pb - P U ) 4 J On 


dp L 
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(Or - Ol) 4 


( 30 ) 

-dtp 

(31) 


In addition, it should be noted that /*f* = Each of the double integrals 
in (23)-(31) is associated with a simple second order polynomial w hose closed 
form evaluation is trivial and is omitted here. 


5 Dyadic Green’s Function 


To complete the evaluation of all matrix elements in ( 19) we must now look 
at B tJ and D lJ% which involve the evaluation of a boundary integral. The 
boundary integral requires numerical evaluation and central to this task is 
an efficient computation of the dyadic Green's function, G e 2 - This Green's 
function introduced in (6) was derived from vector wave functions using the 
procedure espoused by Tai [4]. The resulting modal solution although exact 
is extremely costly to compute, and a more efficient evaluation must be 
presented for practical purposes. This is addressed below. 

Recall from (20) that the boundary integral is of the form 


l B b = j Ili(n. 


<z ).::p,,o l .z l } 


p(a,o,z) x J (7 f2 («-0, z\a,<!>,z') 


y.p(a. o .z ) ■ ii;,(«. o . z : Pj,<pj, ~j) 


do dz dodz 


( 32 ) 


where (a. b) 6 {<?*-} corresponding with the vector component of the test- 
ing and source weights, respectively. Carrying out the vector operations we 
obtain 

if: = “ /, f H';(n, o. Pi, 0 t , z,)\V : (a. o', z' : pj, o,, z s ) 
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(33) 


G<t>o ( a « O — O . r — r )d<p Jr dQclz 


I B! 

1 zo 

= IV. (a. 0 , z\ p t> 0 „ 0 , r'; pj, <j> } , z } ) 

J S ax J S a} 


I 

G<f, z (a , 0 — 0 . z — z )d<p dz d<pdz 
- L L ll '*( rt ’ <?, pi, <j>{, 5,)fF.(a, 0 ', Pj, 0 J, ij ) 

J ^«i J ^aj 

(34) 

1*1 

l o<t> 

G 0 z {a, 0 — 0 — z )clo dz d(j)dz 

= - lFfl>(a, 0 , z;p„ z,)W 0 (a,o .z'-'P^Qj.Zj) 

(35) 


('zz(<t* O — 0 * z — r )d 0 dz d 0 dz 

(36) 

with the weight functions as given by (17). The convolutional nature of 
the Green's function is explicitly shown in (33)-(36). Formally, the Green’s 


function satisfying ( 6 ) and the radiation condition for the cylinder can be 
expressed as 


=) = — — f rm 

v *)\ kJ - AK ) x h^\x) 


(ft . O, z) = — ! — y 1 [ 

( 2 7r ) 2 ^ y_ 

V > 71= — <>D 




Jlnc , - ki :) dk : (: is ) 


G^fo, 0 , c ) = 


— — y r - 



( nk 3 N 

\ 2 tf< 2 , u-)l 

_H ( n\x) 

\ ^ O j 

' Hn {2) (. r)_ 


0 j{no-k z z 


lS) dk z { 39) 


with 0 = 0 - 0 ',= = ;- z\ k p = Jk* - *?, x = fc p a, //‘ 2 >(-) is the n (/l -order 
Hankel function and // n * 2) (-) denotes the derivative of the Hankel function 
with respect to its argument. It is apparent that the evaluation of G' : , G 01 
and G'° using the expressions in (37)-(39) is extremely expensive due to 
the slovs convergence ol both the series and t lie integral especially for near 
self-cell (©, ~ % 0 ) evaluation. 

Bird [7] suggested a more efficient evaluation of the modal Green's func- 
tion by introducing more rapidly convergent integrals through a procedure 
attributed to Duncan [ 8 ]. In doing so, he exploited the symmetry of the 
infinite series, and the indefinite integral was converted into two other inte- 
grals. One of these is definite and the other is a rapidly converging indefinite 


N 




integral. His accelerated modal representations are 


G :: (a , o, ~) = 


= a £ t,«cos(»o) [ Gi[z,n,t)dt+j [ G'i(:.:.- 

f 6 2 (c, n,t)dt + j f G 2 , 

Jo Jo 


G' 00 (n, <5, r) = — ^ e " cos(no) 


6' 0 ‘'(a, <6, =) = 52 s in(nO) f G 3 (z,n,t)dt + j [ G 3 {~.:.~ 

~ a ^o L y ° y o 


j<)^ 

j t)dt\ 

jt)dt 


where 


and 


(40) 

(41) 

(42) 


G](z*nJ) — 

(T/2 0 — 


- j k,jZt 


\i(kaVl-t 2 ) 


G 3 (z, n. t) - 


1 -' k ° st 

1 

(*«) 2 (1 -t 2 ) 

U/2(Wl -t 2 ) 

( nt ) 

1 

\ A cr \/l — t 2 ) 

A' n 2 (Wl - t 2 ) 

ntt~ jk ° st 

1 


+ 


( A'rt ) 2 ( 1 — t 2 ) N^ihas/r^J 2 ) 


(43) 


(44) 

(45) 


-'/ n 2 (o = ^(o + >; 2 (o 

v 2 (/) = -/,?(/) + v;; 2 «) 

€ u = 2 n > 0 

= 1 n = 0 


(46) 

(47) 

(48) 


in which J n (t) and ),,(<) are the usual n 4/, -order Bessel functions. Further 
details regarding refinements required to manage the singularities in the def- 
inite integrals are given in [7], 

Although the accelerated modal solution ( 40 )- ( 42 ) given above are more 
efficient than the original modal form (37)-(39). for large cylinders with ob- 
servation and source points far from each other it is instructive to revert to 
some asymptotic evaluation of ( 37 )- ( 39 ) . Several research groups in the past 
have developed such evaluations of ( 37 )- ( 39 ) , and the most notable of these 
are attributed to Pathak [9], Boersrna and Lee [10], and Bird [11]. These 
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asymptotic formulas are based on similar derivations and differ only in the 
level of approximations offered by each expression. Since the formulas pre- 
sented by Bird are generally the most accurate, we shall concentrate on them. 
However, we note that Boersma and Lee's expressions are particularly well 
suited for sell-cell evaluations and Bird uses part of their solution to improve 
his formula as the source and observation points approach each other [11]. 

The asymptotic evaluations are only useful for large radius cylinders since 
they employ Watson’s transformation [12] to convert the series in (37)-(39) 
into a contour integral. After some manipulation, the contour integral is 
evaluated as a residue series. Each residue contribution is associated with a 
creeping wave which encircles the cylinder from the source to the observation 
while traveling on the cylinder's surface along a geodesic path. Typically, for 
large radius cylinders the contributions of those creeping waves which encircle 
the cylinder once or more are neglected since they are weak in comparison to 
the direct creeping wave contribution. Bird [11] uses a uniform asymptotic 
expansion of the Hankel functions in (37)-(39) and a steepest decent path 
evaluation of the Fourier integral to achieve his formula. He found that 

r) — ~~^ ( l e ~ jk ° S {cos z Q + q{ \ — q ){2 — 3co$ 2 0)) i>( 3 ) 

+ "[ (fT"'" ~ id + (s5“ rAJ) 

( — sm'O + — ) Vjiii) + {.!<■, n ..1t 

-<l 2 d (y + | tan 2 0 - cos 2 0 j v'(i3) (49) 

G° s (a, <?, 5) ~ J -h qe -J k °* s i n 0 cos Q ( i _ :ty( i _ q )) r (;j) 

+, (;(-<•(- -ij)rW)+(ii- §«€»*) J) 

( — j v 2 {3) + ^Csec 2 6>t>oi(.i) 

> f / 1ST 5 ., \ , . 1 

+<r J (Yf ~ ^ sec °) v (/*) (so) 
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G^la, 0 . r) ~ — ^—qt jk<jS (^sin 2 9 + q(l — q)(2 — 'Isin' 2 !))^ v(3) 

+<l (("(•*! - v{3))sec 2 0 + {^tan 2 0 - r(J) 

+ (^g 52 '” 2 ^ ~ ^ an ^) Ul (/^) 

+ ( jt twri^Q - — sin 2 0^j v 2 {3) + ^/3tan 2 0v ol (3) 

+q\3 j^nec 2 0u' (3) + 

w ~ I 2 ? ° + HT Sm °) v (J,) JJ (51) 

The geodesic path length is given by 

trajectory is 0 = tan -1 which is 

shown in figure 7. Also, $ = o or 4> = 2ir — <p depending on which of the 
two diiect paths are taken as illustrated in figure i. In (49)-(51), u(3) and 
t'n{3) represent the* soft and hard surface Fock functions, respect] eel v. These 
functions are characteristic of the creeping waves on a circular cylinder and 
are discussed in detail by Logan [13]. The steps involved in the derivation of 
(49 )— ( 51) are described in [11]. 

Figure 8 illustrates a comparison between the asymptotic formula (49)- 
(51) and the accelerated modal solution (40)-(42). I he given curves corre- 
spond to the magnitude of the dyadic components along a 0 = 10° trajectory 
on an a = 4A cylinder and the dynamic range indicated confirms our asser- 
tion that geodesic paths which encircle the cylinder one or more times need 
not be considered since the resulting large path length results in negligible 
contribution. Figure 9 shows the phase error between the modal solution and 
the asymptotic formula. As seen, the asymptotic formula has small error in 
both magnitude and phase even when observation points are quite close. In- 
deed. Bird has used these formula to investigate the mutual admittance of 
an aperture on a cylinder with minimal difference from a more exact evalu- 
ation [14], It goes without saying that since the computation of the mutual 
admittance involves integration through the singularity of the Green's func- 
tion, we may expect excellent results on using the given asymptotic formula. 
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where q 


and i = ks 


/c 2 + («< I>) 2 and the geodesic 




However, we may find it advantageous to use the modal solution for self-cell 
contributions and the asymptotic formula for all other situations. 


6 Future Tasks 


We have presented a FEM-BI formulation appropriate for aperture anten- 
nas on a metallic circular cylinder. The FEM-BI equation was derived from 
the vector wave equation and converted into a system of equations using 
Galerkin’s procedure with vector-valued weight functions. These functions 
comprise the edge-based elements which are derived from node-based ele- 
ments using Whitney s formalism. Since these elements have a high degree 
of geometrical fidelity for cylindrical arrays and are divergenceless, they are 
well suited for our purposes. We have shown that the boundary integral may 
be numerically evaluated in an efficient manner using asymptotic formulas 
for the dyadic Green’s function while reserving the possibility of using the 
exact, accelerated modal Green s function for the self-cell, if necessary. 

Future work will entail the implementation of the proposed FEM-BI for- 
mulation given in this report. We will exploit the convolutional nature of the 
boundary integral (33)-(36) in the context of the Conjugate Gradient-Fast 
Fourier Transform (CG-FFT) solution technique to maintain low computa- 
tion and 0(N) memory requirements. The FEM portion of the matrix is of 
comse sparse due to the local nature of the finite element method. Upon vali- 
dation of the implementation by comparison with limiting cases and possibly 
measured data, we will undertake a thorough investigation of the properties 
of conformal arrays on a circular cylinder including the following studies: 
mutual impedance, pattern synthesis, element shadowing, and scattering re- 
duction techniques. 

The txpeiience gained in developing an accurate model ol conformal ar- 


rays on a circular cylinder will allow us to extend the analysis to doubly 
cmved cylinders. This will be achieved by employing Geometrical Theory of 
Diffraction (GTD) principles to the Green’s function presented in this report 
so that is is accurate for convex cylinders. In addition, a major challenge will 
be the reduction of memory requirements for the {D tJ \ submatrix (20) since 
for non-circular cylinders the boundary integral will no longer be convolu- 
tional and thus the matrix will not be Toeplitz. We are interested in using 
wavelet transformation techniques to convert the fully populated matrix into 


IS 




an equivalent sparse matrix. An efficient and accurate technique for mod- 
eling conformal arrays on doubly curved surfaces would be very valuable to 
the antenna design community and will allow full utilization of this versatile 
antenna. 
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Figure 8. Magnitude of the asymptotic formula and the modal 

Green s function for a 4A cylinder along a 0 — 10° trajectory. 








